Droplet and cluster formation in freely falling granular streams 
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Particle beams are important tools for probing atomic and molecular interactions. Here we demon- 
strate that particle beams also offer a unique opportunity to investigate interactions in macroscopic 
systems, such as granular media. Motivated by recent experiments on streams of grains that ex- 
hibit liquid-like breakup into droplets, we use molecular dynamics simulations to investigate the 
evolution of a dense stream of macroscopic spheres accelerating out of an opening at the bottom 
of a reservoir. We show how nanoscale details associated with energy dissipation during collisions 
modify the stream's macroscopic behavior. We find that inelastic collisions collimate the stream, 
while the presence of short-range attractive interactions drives structure formation. Parameterizing 
the collision dynamics by the coefficient of restitution (i.e., the ratio of relative velocities before 
and after impact) and the strength of the cohesive interaction, we map out a spectrum of behaviors 
that ranges from gas-like jets in which all grains drift apart to liquid-like streams that break into 
large droplets containing hundreds of grains. We also find a new, intermediate regime in which 
small aggregates form by capture from the gas phase, similar to what can be observed in molecu- 
lar beams. Our results show that nearly all aspects of stream behavior are closely related to the 
velocity gradient associated with vertical free fall. Led by this observation, we propose a simple 
energy balance model to explain the droplet formation process. The qualitative as well as many 
quantitative features of the simulations and the model compare well with available experimental 
data and provide a first quantitative measure of the role of attractions in freely cooling granular 
streams. 

PACS numbers: 45.70.-n, 47.60.Kz, 37.20. +j, 51.10.+y 



Granular material consists of macroscopic, solid parti- 
cles interacting predominantly via contact forces p~H3] . In 
modeling granular flows these interactions are typically 
taken as purely repulsive. However, there are important 
circumstances, from agglomeration in fluidized particle 
beds to dust accretion in proto-planetary discs, where 
attractions can compete with the particle weight or with 
forces produced by particle collisions and then lead to 
the formation of stable granular clusters [4-7 . In dry, 
nominally free-flowing granular material the attractions 
are small and short-ranged, and are usually associated 
with van der Waals forces or capillary bridges due to 
a few layers of adsorbed molecules [TJ |S]. To measure 
the resulting residual cohesion between grains, atomic 
force microscopy (AFM) has proven useful [4, 9, 10 , but 
this technique mimics the static limit of central, head- 
on collisions. Under dynamic conditions, energy is also 
dissipated due to the inelastic nature of collisions, includ- 
ing rolling and sliding friction during non-central impacts 
[5l HTHT4] . Simulations have started to address the com- 
petition between inelasticity and cohesion during colli- 
sions in freely cooling granular systems, but so far have 
focused on the limit of a dilute granular 'gas' [6} fT3| fT5] . 
Little is known about the complex dynamics that lead to 
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clustering in the dense limit, where many weakly cohesive 
particles collide and interact in rapid succession. 

New possibilities to tackle this problem have emerged 
from recent experiments on freely falling granular 
streams [T6HT8] . In these systems, a dense stream of 
particles emerges from a small opening at the bottom 
of a reservoir and is accelerated downward by gravity. 
During vertical free fall, particle interactions lead to spa- 
tial inhomogeneities that can be detected downstream. 
This makes it possible to observe the effects of attrac- 
tive forces as small as nanoNewtons between macroscopic 
grains [18] • Such forces constitute an effective 'surface 
tension' that can affect the bulk dynamics of the stream. 
In principle, the sensitivity afforded by such experiments 
should provide a means to delineate cohesive from purely 
collisional contributions to structure formation. This is 
because in the accelerated, co-moving frame the stream 
undergoes extensional flow and particles initially touch- 
ing will eventually separate along the axial direction un- 
less they are held together by attractive forces. This is 
similar to what happens during free expansion of molec- 
ular jets p~9j|20]. The details of cluster evolution, shape, 
and size therefore contain information about the inter- 
particle interactions. However, to access this informa- 
tion, a connection between cluster properties and rele- 
vant particle parameters must first be established. 

As a first step, we report here on a set of systematic 
simulations designed to differentiate between the roles of 
inelasticity and cohesion in driving droplet formation or 
clustering. Such simulations can give access to aspects 
difficult to probe experimentally, such as the local par- 
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ticle configurations and dynamics in the interior of the 
stream. Furthermore, they can explore parameter ranges 
that so far have not been tested in experiments. 

Covering large regions in the cohesion-inelasticity pa- 
rameter space, we not only reproduce qualitatively the 
observed experimental findings, but gain new insights re- 
garding the roles of inelasticity and cohesion as well as 
the mechanism through which droplet formation occurs. 
In terms of the gross behavior, we are able to delineate 
regimes ranging from sprays of isolated particles to ag- 
gregation into clusters by collision and capture to break- 
up into droplets as in Fig. [TJd) or even larger, solid-like 
chunks. On a much more local scale, by tracking the aver- 
age contact number of particles in the stream, we demon- 
strate that clustering and droplet formation are signa- 
tures of attractive forces, while inelasticity plays only a 
secondary role by helping to collimate the stream. We 
also show that the average network size correlates directly 
with the cohesion strength. This opens up unique possi- 
bilities to quantify cohesion by simply measuring the net- 
work size and, more generally, to use granular streams as 
sensitive probes of inter-particle attractive interactions. 



I. RESULTS AND DISCUSSION 

We model the interactions between particles with 
molecular dynamics, including all three-dimensional de- 
grees of freedom, following the approach described in 
Refs. [23H25] , The system starts as an ensemble of up to 
70,000 monodisperse spheres of diameter d inside a cylin- 
drical reservoir (hopper). Pulled by gravity, grains leave 
the hopper bottom through an aperture of diameter Do? 
and the hopper is continually replenished at the top to 
keep the average fill height constant. For practical pur- 
poses, we work primarily with d = 200 jam spheres and 
aperture size Do = 3 mm (Figs.[2][9|, although we present 
data with d = 100 /im and Do = 2 mm in Fig. [T] Com- 
pared to the simulations in Fig. [TJ the smaller stream di- 
ameter to particle diameter ratio D$/d = 15 for Figs.|2][9] 
is computationally more efficient, but, as in the experi- 
ments [18], it leads to more scatter and less definition in 
the resulting stream features. Since the larger grain mass 
increases the collisional kinetic energy, the strength of co- 
hesion required for droplet formation is also increased. 

Normal and tangential contact forces are modeled with 
a linear spring-dashpot force law. Collisional and fric- 
tional dissipation are parameterized by coefficients of 
restitution e and static friction \i. Cohesion between 
grains is incorporated as a constant attractive force F co h 
which acts over a distance l c « d. This force only turns 
on after particle surfaces have come into contact. As a re- 
sult, completely separating two particles initially in con- 
tact costs an amount of energy W co h = J^ lc Fcohdr = 
F co hl c . This work has the same origin as W a d introduced 
by Brilliantov et al. [TT] in their model for viscoelastic 
adhesive collisions, although in our case we restrict our 
particles to a bare coefficient of restitution that is ve- 
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FIG. 1: (Color online) Snapshots from simulations of granular 
streams consisting of d = 100 /am spheres freely falling out 
of a Do = 2.0 mm aperture, (a)-(d) Coefficient of restitution 
e = 0.40 and cohesive strength F coh = 100 nN. (e)-(h) e = 0.88 
and F co h — 1 nN. The images follow the grains in the co- 
moving frame and are taken just below the aperture [(a),(e)], 
and at distances z = 0.7 cm [(b), (f)], 1.9 cm [(c), (g)] and 
z — 3.6 cm [(d), (h)] from the aperture to the top of each 
frame. Grains are color coded according to their rms velocity 
difference with neighboring grains within a distance of 1.5gL 
Isolated grains are shown in pink. Images generated via |21j . 
For movies see |22j . 



locity independent and a constant cohesive force. Physi- 
cally, this model is often used to mimic short range forces 
such as capillary bridges and van der Waals interactions. 
For a thorough discussion of these forces and their con- 
sequences, we refer the reader to the review by S. Her- 
mingaus [26] . 

We do not include interactions between the particles 
and an interstitial fluid (air) since experiments performed 
under different levels of vacuum showed that the basic 
features of droplet formation do not depend on the pres- 
ence of air drag [18]. For the results reported here, we 
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set p = 0.5, l c = 100 nm and the specific grain density to 
that of glass (p = 2.5 kg/m 3 ), and we explore the effects 
of collisions and cohesion by varying e and F co h- 

The results shown in Fig. [T] for d = 100 pm grains 
falling from & D Q = 2 mm aperture reproduce the two 
types of behavior observed in experiments for similarly 
sized glass spheres with short-ranged cohesive interac- 
tions of comparable strength (as was confirmed by AFM 
[18]), namely spraying and droplet formation. For a 
highly cohesive, highly dissipative stream, as shown in 
Fig. [TJa)-(d), the stream begins by narrowing in width 
as is commonly observed in liquid streams. Further 
downstream undulations form and eventually evolve into 
well-defined droplets connected by necks. These necks 
thin and eventually break. When the cohesive strength 
was reduced in the experiments by roughening the parti- 
cle surfaces (also confirmed by AFM), droplet formation 
ceased and the stream behaved more like a collimated 
spray. The same qualitative change is reproduced when 
the cohesive strength is decreased in the simulations [ 
Fig. [TJe)-(h)]. Fig. [TJa)-(d) also shows that droplet for- 
mation goes hand in hand with the development of spatial 
heterogeneities in the relative particle velocities: regions 
of larger relative velocity deform into necks and stay 'hot' 
even for a short time after they break, while the interior 
of the droplets is 'colder' and they appear to be essen- 
tially frozen once the separation process is complete [for 
details of this evolution see supplemental movies [22]]. 

The overall shape of the droplets can be characterized 
by their length to width aspect ratio Xd/wd- For the 
purpose of establishing shapes, height and width data 
was gathered for all droplets of 30 particles or more in 
simulations whose cohesive force F co h was greater than 
100 nN (see Appendix for details). As shown in Figs. [2^ a) 
and (c), most droplets have Xd/wd between about one 
and three, in agreement with experiments [T8] , One fea- 
ture the simulations produced that was not seen in ex- 
periments is the tail of droplets whose aspect ratios are 
larger than three [Fig. J2jc)]. These data are most likely 
the result of the inclusion of neighboring droplets that 
did not fully separate before the simulation was termi- 
nated (an observation supported by the slight increase 
at Xd/wd ~ 4). The histogram in Fig. [2|c) shows that 
~ 75% of droplets have an aspect ratio between one and 
three and that the most probable aspect ratio lies very 
near to two. As can be seen in Fig. [2^b) , for a given in- 
elasticity the length and width of droplets increases with 
the amount of cohesion. 

These observations are in contrast with inviscid liquid 
streams undergoing Rayleigh-Plateau breakup. First, liq- 
uid breakup has unstable wavelengths only for Xi/Dq > ir 
(here we use Xi to distinguish the liquid lengthscale from 
the granular lengthscale) [27]. It is important to note 
that the definitions of these aspect ratios are slightly dif- 
ferent, with Dq (the local unperturbed jet diameter) as 
the denominator in the liquid case and Wd (the droplet 
width) in the granular case. However, if we make the 
more direct comparison and calculate Xd/Do for the 



granular case, we find that Xd/Do ~ 2/3, drawing an 
even clearer distinction between the granular and liquid 
streams. This comparison is still complicated by the fact 
that the granular stream is continually being stretched 
axially by gravity, and if we account for this by not- 
ing that the length of fully formed droplets is necessarily 
longer than the corresponding lengthscale at the aperture 
(which we will call Ao), we deviate still further from the 
liquid case arriving at Xq/Dq < 2/3. Second, for inviscid 
liquids undergoing Rayleigh-Plateau breakup, the size of 
droplets is independent of the surface tension. In fact, the 
only parameter governing the size of droplets is the aper- 
ture diameter, with Xi > ttDq determining the minimum 
unstable wavelength and Xi « 4. 5 A) yielding the fastest 
growing wavelength [27 . As seen in Fig.^b), the droplet 
size does depend on the cohesive strength in the granular 
case. On a microscopic level, surface tension arises from 
cohesive bonds between constituents, whether they are 
molecules or granular particles, i.e. 7 ~ W co h/d 2 (where 
7 is the surface tension, W co h — F co hl c is the cohesive 
energy stored in a bond, and d is the particle diameter). 
Even so, these observations make it clear that the mecha- 
nism by which this "surface tension" leads to systematic 
breakup in the granular case differs from the liquid case. 

Along with cluster aspect ratio, our simulations also 
reproduce the observed linear growth of the separation 
between adjacent droplets with time [Fig. [3]. Identify- 
ing the grains which make up a fully formed droplet, 
we measure the vertical distance between the centers of 
mass for adjacent droplets Az(t) and compute the sepa- 
ration velocity Av(t) = d(Az)/dt both in the free falling 
stream (t > 0) and inside the hopper (t < 0). Even in- 
side the hopper Az > and Av > 0, indicating that 
despite the large amount of shear there is relatively little 
vertical mixing. As pairs of droplets approach the aper- 
ture they rapidly accelerate apart until they leave the 
aperture with a separation velocity Avo. After leaving 
the aperture the increase in Av abruptly stops, though 
the droplets can still interact with each other, either ac- 
celerating or decelerating the separation of neighboring 
pairs. Contrasting Fig. [3] (a) and (b), it becomes appar- 
ent that droplets in streams with larger cohesive forces 
undergo stronger, more sustained interactions, whereas 
droplets in the streams with the smaller cohesive force 
exhibit weaker, more erratic interactions. These interac- 
tions could not be resolved in [18] since experimentally 
Az could only be measured at later times after well de- 
fined undulations had formed. After a droplet and its 
neighbors are no longer connected to other droplets Av 
remains fixed at a final separation velocity Avf. 

The simulations allow for the systematic investigation 
of parameters which are not easily changed in experi- 
ments, such as F co h and e. Fig. [4] represents a 'phase 
diagram' based upon these two parameters that delin- 
eates different regimes according to general, macroscopic 
types of stream behavior: spraying, clustering, droplet 
forming, and dripping. The images provided for sev- 
eral of the (F co h,e) pairs show snapshots at a depth 
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FIG. 2: (Color online) Droplet shapes for d — 200 fim spheres freely falling from a Do — 3.0 mm aperture, (a) Droplet height 
Xd vs. width Wd for simulations in droplet forming regime. Only a sampling of simulations is shown here to avoid clutter [data 
for all droplet forming simulations are included in Fig. [5Jc)]. Inset demonstrates what is meant by width Wd and height A^. 
(b) Droplet width Wd (red) and height A^ (blue) vs. F co h for e = 0.61. Error bars show standard deviation of measured Ad and 
Wd- Aperture diameter is indicated by dashed line, (c) Distribution of aspect ratios Xd/wd for all droplet forming streams. 



z = 4.0 cm downstream. This depth was chosen to al- 
low sufficient time for some of the salient stream features 
to emerge (c.f. Fig. [I]). The images highlight the differ- 
ent regimes of stream behavior, but also reveal that it 
is possible for streams in different regions of phase space 
to appear nearly indistinguishable by shape alone, such 
as Figs. 4ja)-(c) and Figs. |4^e)-(f). In order to tease 
out subtle differences in the structure and connectivity 
of these streams we compute each grain's contact num- 
ber C, defined as the number of neighbors whose centers 
are less than d + l c away from the grain's center. This 
includes both surface-to-surface contacts and neighbors 
within range of the cohesive force F co h. For example, 
Figs.ga) (F coh = 0, e = 0.88) and gb) (F coh = 100 nN, 
e = 0.88) both appear to be diffuse sprays based on the 
2-dimensional projection of particle positions alone, but 
the color coded images reveal a substantial number of 
contacting pairs for the 100 nN stream and virtually none 
for the zero cohesion stream. We refer to streams that 
evolve into isolated particles as 'sprays' and distinguish 
this behavior from the formation of small aggregates in 
the 'clustering' regime. As an example, even though a 
stream corresponding to e = 0.04 and F co h — [ Fig.gf)] 
appears much denser and more collimated, there are only 
a few scattered contacting pairs. In Fig. ge) (F co h = 50 
nN, e = 0.40), however, most grains have at least one 



contacting neighbor and some have as many as four or 
five. This already indicates that, while the inelasticity 
strongly affects the degree to which a stream remains 
collimated, a non-zero cohesive force is essential for main- 
taining contacts between adjacent grains. Coloring the 
data points in the phase diagram by the average contact 
number (C(z)) at z = 4 cm reveals no sharp boundaries 
between different regimes but rather gradual crossovers. 

Towards larger cohesive forces and smaller coeffi- 
cients of restitution lies the droplet forming regime 
[c.f. Fig. EJh)] discussed already in connection with 
Figs. [TJa)-(d). The behavior here is reminiscent of liquid 
flows from a small opening and individual droplets can 
contain hundreds of grains. If we increase the cohesion 
even more and approach the right-hand edge of Fig. [4j 
we enter a regime in which there is no longer a steady 
outflow from the aperture. In this 'dripping' regime large 
aggregates of grains slowly ooze through the aperture [ 
Fig.gd)]. Once the weight of these aggregates becomes 
too large, big chunks break off and pull away. Further 
increasing the cohesive force, the grains become perma- 
nently stuck inside the hopper and all flow stops. 

Since droplets and clusters are aggregates of grains 
with vanishing relative velocity, it is tempting to use 
the notion of a granular temperature to characterize the 
roles of inelasticity and cohesion. In the kinetic the- 
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FIG. 3: (Color online) Droplet separation dynamics. Data 
are for the same grain diameter and aperture size as in Fig. [2] 
Distance between the center of mass Az{t) (insets) and sep- 
aration velocity Av(t) = d(Av)/dt for four pairs of adjacent 
droplets from streams with e = 0.61 and (a) F co h = 400 nN 
and (b) F co h = 1100 nN. For each pair, the the time t is off- 
set so that at t = the center of the pair of droplets passes 
through the aperture (z — 0). Images to the right show the 
pairs of droplets tracked here at a depth z — 18.4 cm below 
the aperture (t ~ 180 ms). 



ory of granular flows, the granular temperature is ob- 
tained by computing local deviations of grain veloci- 
ties from a coarse-grained mean velocity field V(x, t) 
[2j [28]. In freely falling streams, however, the veloc- 
ity field varies strongly in both time and space. Nev- 
ertheless, to capture the variations in relative veloci- 
ties we can calculate a mean squared velocity deviation 
<(Av) 2 ) = ((A^) 2 ) + ((Avy) 2 ) + ((Av z ) 2 } of each grain 
relative to its neighbors within a small region. Below we 
will loosely refer to this measure as a granular temper- 
ature, but it is important to remember that it contains 
contributions from collisional velocity fluctuations and 
as well as from gradients in the mean velocity across a 
grain's neighborhood. 

In Fig. [5] we plot the mean squared velocity devia- 
tions versus depth for a sampling of parameter values 
from the representative regions of the phase space. In- 
side the hopper and sufficiently far above the aperture, 
velocity deviations are small and grains move in unison. 



As the aperture is approached from above, shear 'heats' 
up the granular medium, leading to a pronounced peak 
in ((Av) 2 ). At this point, the rms velocity deviations 
reach a value that is a significant fraction of the mean 
outflow velocity, (v z ) « 10 cm/s. We find that the max- 
imum deviations are nearly the same for a wide range of 
cohesive forces and inelasticities, suggesting that the dy- 
namics inside the hopper and the initial flow conditions 
right at the aperture are largely unaffected by changing 
these parameters. 

Once the particles leave the hopper the 'cooling' phase 
begins and differences emerge. Figs. 5jb) and (c) show 
separately the transverse and axial contributions to 
((Av) 2 ). For weakly cohesive streams, and irrespective 
of e, the large- z behavior is seen to approach a power 
law for both contributions, albeit with the exponent of 
the transverse component being slightly smaller. Streams 
that undergo strong droplet formation (green and blue 
traces) show deviations from a simple power law, exhibit- 
ing first an excess of fluctuations (c.f. the shear visible 
in the neck regions in Fig. [I]) followed by a rapid drop of 
the local droplet temperature after break-up, as cohesive 
forces bind grains together. 

In the absence of drop formation the observed axial 
cooling follows, to very good approximation, the behav- 
ior expected just from the stretching of streams of non- 
interacting particles, i.e., simply from the velocity gradi- 
ent due to gravity. As particles fall out of the aperture, 
they acquire a free fall velocity 



v z {z) = v z $^l + 2gz/v 2 z ^ 



(1) 



where v z $ is the axial velocity at the aperture and g is the 
gravitational acceleration. The axial velocity deviations 
over a local range Az are then computed as 



((Av z (z)) 2 ) = ((v(z + Az)-v(z) 2 ). 



(2) 



This predicts not only the 1/z decay of the axial tem- 
perature at large z, but, as seen in Fig. [5jd) (dashed 
line), fits the data of non-droplet-forming streams quan- 
titatively from the aperture on down (in these fits the 
initial velocity v z $ is known and the only free parameter 
is the neighborhood size Az). Since z is proportional to 
the square of the falling time, the axial temperature in 
Fig. km d) drops as 1/t 2 , coincidentally the same behav- 
ior as predicted for the early stages of a non-expanding 
granular gas [6]. We note that this differs from the 1/z 2 
decay proposed by Amarouchene et al [16J based on fluc- 
tuation measurements at the surface of streams falling in 
air. 

While the mean squared velocity deviations reveal 
the important role played by gravity in stretching and 
thereby diluting and cooling the stream, they are not 
particularly sensitive to details introduced by inelastic- 
ity or cohesion. We see this directly in Fig. [5] where 
temperature curves for a range of F co h and e values ef- 
fectively collapse on top of each other except when the 
stream undergoes gross changes in the behavior, such as 
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FIG. 4: (Color online) Phase diagram parameterized by cohesive force F co h and restitution coefficient e. Data are for the same 
grain diameter and aperture size as in Fig. [5] Data points are color coded by the average contact number (C(z)), computed at 
depth z = 4 cm. Images (a)-(h) from the simulations at the same depth are shown for selected (F co h, e) pairs. Grains in these 
images are colored by their local contact number according to the same color scale. Data along the left vertical axis corresponds 
to F co h = 0. The image in the upper right highlighted by the red border corresponds to z = instead of 4 cm. 



the onset of large droplets. It is therefore important to 
explore additional and possibly more useful indicators of 
the roles of F co h and e. 

One of these indicators is the depth dependence of the 
average contact number, (C(z)}. To emphasize changes 
relative to the initial particle configuration near the aper- 
ture, we normalize (C(z)) by the average initial contact 
number, (Co), measured at z = 1 mm. For cohesive 
streams in the clustering, droplet forming, or dripping 
regimes, the normalized contact number always increases 
with z [ Figj6ja)]. Contrarily, in non-cohesive streams, 
(C(z))/(Co) always decreases, independent of the inelas- 
ticity and even in streams with restitution coefficients 
approaching zero. Since the interparticle distance, and 
hence the collision rate, in the stream steadily decreases 
with depth, the usual pressure gradients associated with 
inelastic collisions are rendered less effective, and there is 
no restoring force to keep contacting grains together. As 
a result, even the slightest separation velocity-essentially 
ensured by the velocity gradient-will eventually cause a 
cluster of contacting grains to drift apart. Only in co- 
hesive streams is it possible to form particle aggregates 
that will remain stable against small velocity differences 
in the absence of further collisions. 



While the normalized average contact number draws 
a clear distinction between cohesive and non-cohesive 
streams, it does not differentiate cleanly between the 
more subtle consequences of cohesion. Such differences 
are highlighted when we go beyond nearest neighbor con- 
tacts and look at the full network of particles connected 
to a given grain. We do this by defining the network 
size (N(z)} as the average number of grains connected 
to a grain at a distance z below the aperture (see Ap- 
pendix for averaging procedure). We define two grains 
as connected if a path exists between them that consists 
of pairs of grains whose center-to-center distances are no 
larger than l c + d. Fig. [6^b) shows the results for streams 
with e = 0.61 and a range of F co h values. 

The behavior of (N(z)} can be classified into three ba- 
sic categories. For F co h = we find that (N(z)} mono- 
tonically decreases toward unity, in agreement with the 
conclusions drawn from the evolution of (C(z)). This is 
the spraying regime. For F co h less than about ~ 100 nN, 
(N(z)} first dips sharply as particles leave the hopper, 
but then undergoes a prolonged period of slow growth, 
after which it saturates at values less than ~ 30 grains. 
We interpret this as the signature of clusters that form by 
collide-and-capture events. The dynamics in this regime 
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FIG. 5: (Color online) Average local velocity deviations as a function of depth. Data are for the same grain diameter and 
aperture size as in Fig. [2] The average axial velocity at the aperture was v z ,o « 10 cm/s. (a) Mean square velocity deviations 
((Av) 2 ) inside the hopper [z < 0) and in the falling stream [z > 0). (b) Radial velocity deviations ((Av r ) 2 ) = ((Av x ) 2 + (Av y ) 2 ) . 
(c) Axial velocity deviations ((Av z ) 2 ). Data in (a)-(c) correspond to: F co h = 0, e = 0.88 (orange), F co h = 0, e = 0.48 (pink), 
F coh = 50 nN,e = 0.40 (purple), F coh = 100 nN, e = 0.88 (red), F coh = 400 nN, e = 1.0 (black), F coh = 300 nN, e = 0.61 
(green), F co /i = 1000 nN, e = 0.4 (blue), (d) Same as in (c) with droplet forming streams (green and blue curves) excluded. 
Dashed line: fit to Av 2 (Eq. 2), giving Az = 200 /mi. 



consist of an intricate interplay between the effects of the 
velocity profile at the aperture, the stretching and dilut- 
ing of the stream via gravity, and low energy collisions 
which permit initially separated grains to become perma- 
nently connected. The capture process can start once the 
stream has fallen and cooled to a low enough level such 
that the cohesive energy F co hl c exceeds the pre-collisional 
kinetic energy [5j |6] . This does not continue indefinitely, 
however, because the stretching due to gravity drives the 
collision frequency toward zero. This situation is remark- 
ably similar to the trade-off between cooling and dilution 
observed in free expansion experiments on cryogenic He 
jets [19]. 

For larger F co h, between 100 nN and 2000 nN, the be- 
havior reverses. Now (N(z)) quickly increases to a peak 
value before settling down to its asymptotic value at large 
z. This is the droplet regime. The initial increase in 
(N(z)} is due to the liquid-like narrowing of the stream 
below the aperture, where the associated radial influx 
and sufficiently large cohesive energy work together to fa- 
vor network growth. The subsequent decrease signals the 
break-up of the stream into droplets. The initial network 



growth for the strongly cohesive streams becomes smaller 
as the cohesive force is increased. At the same time, 
the peak network size saturates at some maximum value 
(~ 1100 grains). We also see that the break-up is more 
pronounced in more cohesive streams, i.e. the sudden 
drop in (N(z)) occurs more rapidly [inset to Fig. |6^b)]. 
Further inspection of Fig. [6jb) reveals that the location 
of this dropoff moves up toward the aperture as F co h is 
increased, suggesting a natural progression from droplet 
forming to dripping streams. 

The qualitatively different shapes of the (N(z)) curves 
reveal a fundamental difference between clustering and 
droplet forming streams. A clustering stream starts as 
'gas' of individual grains or collections of small network 
fragments which grow larger by capturing additional par- 
ticles during collisions. In contrast, a droplet forming 
stream starts as large network of grains that eventually 
breaks into smaller portions, again similar to the forma- 
tion of drops in liquid He jets [20 . A potential crossover 
between these two modes can be identified by tracing 
the asymptotic, large- z limit of (N(z)} as a function of 
cohesive strength. In Fig. [6^c) we plot the saturation 
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FIG. 6: (Color online) Evolution of contact number and network size. Data are for the same grain diameter and aperture size 
as in Fig. [5] (a) Normalized average contact number vs. distance below aperture. Data shown for F co h = nN and e = 0.04 
(solid blue); 0.4 (solid red), 0.88 (solid green) and for F coh = 100 nN, e = 0.88 (blue dots); 300 nN, 0.61 (green dots); 1000 nN, 
0.40 (red dots). Inset: Initial contact number vs. restitution coefficient for F co h = (triangles), 100 nN (squares), 500 nN 
(circles), and 1000 nN (diamonds), (b) Average network size vs. distance below aperture for e = 0.61 streams. From bottom to 
top along right-hand edge: F coh = nN, 30 nN, 50 nN , 70 nN, 100 nN, 150 nN, 200 nN, 600 nN, 1000 nN, 1300 nN, 2000 nN. 
(c) Large- z saturation value of average network size vs. cohesive force for e = 0.61 streams. 



network size (N 3a t), defined as (N(z = 20 cm )) vs. 
F co h. Clustering streams undergo accelerated growth in 
(N sa t) as F co h is increased, (e.g. before about ~ 100 nN) 
while droplet forming streams exhibit a slow turnover in 
which the growth rate decreases. If we inspect Fig. j6^b) 
more closely, we see that the curves undergo the qualita- 
tive change from 'slow growth' (clustering) to 'eventual 
breakup' (droplet formation) around ~ 100 nN as well. 

In Fig. [T^a) we examine how F co h affects the dynamics 
of the stream breakup in the droplet forming regime. The 
average separation velocity between adjacent droplets ex- 
hibits a weak dependance on F co h- The averages of the 
initial (Avq) and final (Avf) separation velocities are es- 
sentially identical, though the distribution of Avf is con- 
siderably wider due to the interactions between droplets 
as they separate. When A^o is plotted against the mean 
number of grains in individual pairs, N pairi the data for 
different F co h collapse onto a single curve [Fig. (7^b)], in- 
dicating that the dependence of (Avq) on F co h stems en- 
tirely from the increased network size at larger cohesive 
forces. The data similarly collapse when plotted against 
the center of mass distance between the clusters at the 
aperture Az = Az(t = 0) [Fig. ^c)]. From Eq. [I] the 
velocity difference of two points in the stream near the 
aperture and separated by a distance Azq <C v^ Q /g is 

A,^f^). (3) 

The fit shown in Fig. [7] yields gd/v z $ = 1.4 cm/s, rea- 
sonably close to the value 1.7 cm/s predicted from just 
the free-fall velocity profile. This shows that the droplet 



separation is primarily driven by the free-fall velocity gra- 
dient along the stream and is set by the initial droplet 
separation, which is on the order of a few grain diam- 
eters. This scale is in agreement with the extrapolated 
value of the length scale measured in previous work on 
granular streams [T7] . 

So far, the simulations have allowed us to reproduce 
experimental results and access data that would other- 
wise be impossible to measure. We can also use them to 
test possible mechanisms for the regions in phase space 
outlined in Fig. [4j For the capture phase, it is natural to 
simply think of the stream as a gas of colliding grains. 
In this scenario, structure formation takes place mainly 
through collide- and-capture events. The picture is sim- 
plified if we assume all collisions are binary. As pointed in 
[TT] , the minimum relative pre-collision velocity v ca p two 
particles must have in order to separate after a collision 
is 

lie F co h / a\ 

Vcap = \ ( 4 ) 

If details regarding the density and the distribution of 
relative particle velocities at the aperture are fairly con- 
stant, it is natural to expect that the mean contact num- 
ber C will depend primarily on v capi i.e. F co ^, Z c , and 
e. This is tested in Fig. [8j Indeed, plotting (C) vs. the 
capture velocity as opposed to the cohesive force does 
bring the data together appreciably even for higher force 
streams, but it does not quite capture the whole story. 
Streams with lower e systematically deviate to higher 
mean coordination numbers. This observation is consis- 
tent with the fact that the introduction of a cohesive force 
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FIG. 7: (Color online) Initial and final separation velocity varying F co h. Data are for the same grain diameter and aperture 
size as in Fig. [2] (a) Average initial separation velocity Avo (black triangles) and final separation velocity Avf (red squares) 
vs F co h. (b) Avo for pairs of adjacent droplets from streams with e = 0.61 and F co h from 400 nN to 2000 nN plotted against 
the mean number of grains in each pair N pa i r = (Ni + Nj)/2 for droplets i and j. (c) Avo vs droplet separation at the aperture 
Azo = Az(t = 0) scaled by the grain diameter d. Solid line is a fit to Avq = a(Azo/d) with gd/v z ,o = 1.4 cm/s. 



leads to an effective coefficient of restitution [TT], in our 
case given by 

It follows that even if two dense granular flows consist 
of particles with the same capture velocity, the one with 
the higher e (and consequently e e ff) should exhibit more 
structure formation. 

Deep in the droplet forming regime, it is tempting to 
look for a lengthscale that falls out from our free param- 
eters since the size of droplets scales with the cohesive 
energy. We can begin to picture a mechanism for the se- 
lection of this lengthscale by recalling two observations. 
The first is that well into the droplet regime virtually 
all grains leaving the aperture are connected in one large 
network. The second is that the velocity gradient at the 
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FIG. 8: (Color online) Contact number vs. F co h and v cap - 
Data are for the same grain diameter and aperture size as in 
Fig. [5] (a) Mean contact number C measured 4 cm below the 
aperture plotted against the cohesive force F co h for e = 1.0 
(green diamonds), 0.88 (blue triangles), 0.73 (orange squares), 
0.61 (red circles) and 0.40 (magenta squares), (b) Mean con- 
tact number at 4 cm from the aperture plotted against the 
capture velocity v cap with same symbols as in (a). 



aperture plays an important role. Putting these together, 
one might suspect that droplet formation is closely re- 
lated to the kinetic energy imprinted on the stream by 
the gradient. As a simplified model, we assume that our 
stream is a well-connected column of grains of infinite 
length and diameter Do, isolated from all external forces 
and subject to an initial axial velocity gradient of con- 
stant value 7o- We partition this stream into segments 
of length Ao which will eventually become droplets. For 
an individual droplet, conservation of energy requires 

K i + W-U i = K f + AE int , (6) 

where Ki is the droplet's initial kinetic energy, W is the 
work done on it by neighboring droplets, U\ is the initial 
energy stored in bonds (Ui > 0), Kf is the final kinetic 
energy, and AE int is the change in the internal energy of 
the droplet. We can calculate the initial kinetic energy 
in the rest frame of the droplet from our knowledge of 
the initial velocity gradient 

Ki = [ X ° /2 \p4>oir(^)\joz) 2 dz = e\l (7) 

J Ao /2 \ / 

where p is the specific density of the grains and 0o is 
the initial packing fraction. In this model with identi- 
cal droplets, the forces on a droplet from the right and 
left must be equal throughout the entire process, so its 
center of mass velocity must remain constant. Though 
not strictly the case in the simulations where the droplet 
size can vary, on average this is true, as evident from 
Fig. [J a). We assume that Kf is identically zero, i.e. no 
energy remains in the form of particle vibrations or bulk 
rotations. 

The work W is more difficult to estimate because it 
depends strongly on the dynamics of the droplet forma- 
tion process. We know it is present and non-negligible 
by the separation dynamics shown in Fig. [3J which also 
suggest that this work grows with F co h. We can estimate 
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its magnitude by noting that it is just the total force in 
the cross section between the droplets, F cs (i), integrated 
over the effective length which the force pulls. In the rest 
frame of our droplet, this work is approximately 



W 



dl 

F cs (t)—dt = 70 A 



£=0 



dt 



t=t b 



F cs (t)dt a\ 



t=o 



(8) 

where a depends on the time to breakup £5, the ve- 
locity gradient 70, and the cohesive force F co h- The 
assumption of a constant velocity gradient is not per- 
fect, since interactions with neighboring droplets and in- 
ternal dissipation can, in principle, alter it. However, 
we have seen from Fig. 5^d) that even droplet form- 
ing streams initially exhibit the same 1/z decay in the 
mean square velocity deviations seen in non-droplet form- 
ing streams, indicating that deviations from a constant 
gradient are only significant close to break up. The 
integral / 



t=o 



F cs (t)dt depends on the full evolu- 
tion of the bond configurations between the neighbor- 
ing droplets, making it difficult to evaluate. However, 
from Figs. [2] and [3] we see that the droplet width, length 
and break up time vary only weakly with F co h-> suggest- 
ing that these changes are minor and to lowest order 
/ oc F co h. In Fig- we demonstrated that the ini- 
tial velocity gradient 70 is independent of F co h- Putting 
these observations together suggests that, to lowest or- 
der, a cx F co h- 

The last quantities we must account for are Ui and 
AE int . These are closely related, and to account for 
them we start with a simple example. Imagine a single 
droplet (with no neighbors) of length Ao and initial ax- 
ial velocity gradient 70. The particles in this droplet are 
perfectly elastic and frictionless. Assuming the total co- 
hesive energy remains constant, the droplet can remain 
intact if the energy stored in bonds exceeds the initial 
kinetic energy, i.e. Ui = N^Wcoh > K% = eAjj, where 
Nb is the total number of bonds in the drop. The num- 
ber of bonds Nb should grow linearly with the droplet 
length (7V5 = n&Ao), which implies that the maximum 
length our droplet can be in order to remain intact is 
Ao = y^n^Wcoh/t' Since we assumed no bonds are bro- 
ken, energy is conserved and the particles will exchange 
kinetic for potential energy as bonds are 'stretched' and 
retracted. If we now assume that during this process the 
number of bonds actually grows [as it does in Fig. [6ja)], 
our stream is able to counteract even more kinetic en- 
ergy (note that creating new bonds adds no kinetic en- 
ergy since our force is hysteretic). We also expect that 
the number of bonds gained should scale linearly with 
the leng th. Thus our dro plet length can now grow to 
Ao = \J (n® + nbg)W C oh/£, where n^g is the number of 
bonds gained per unit length. As more bonds are gained, 
the droplet finds itself in a deeper effective potential well 
and must have more kinetic energy-i.e. more length-to 
escape. The same is true if we add dissipation in the 
form of friction or inelasticity. For each dissipative term 
added, more kinetic energy can be extinguished and this 



produces longer droplets. The work done by neighboring 
droplets has the opposite effect, reducing the maximum 
stable Ao- However, for droplet formation we must al- 
ways have Ui + AEi nt > W. The local rearrangements 
and deformations that lead to this dissipative term are 
easily seen in our supplemental movie [22] , where we show 
the evolution of a single droplet. A similar phenomenon 
occurs in liquids, where increasing the viscosity leads to 
larger droplet masses in the slow dripping from a tube 

Assuming the contact distribution and collision dy- 
namics are uniform across the length of the droplet, dis- 
sipation from friction and inelasticity should also scale 
with the length of the droplet, so that we can write the 
sum of Ui and AE int as 



Ui + AE iT 



(9) 



where (3 includes contributions from all initial and gained 
bonds, friction and inelasticity. Putting our proposed 
forms for W, and Ui + AE int into Eq. [6j we arrive at 
an expression for the initial length of droplets: 



Ao 



P-a 



1/2 



(10) 



Although precisely predicting Ao requires analytical ex- 
pressions for a and j3 as a function of all our simulation 
parameters and initial conditions, we can make some 
progress by just looking at the scaling with F co h- We 
have already argued that a (and consequently W) and 
Ui scale linearly with F co h- Since there is no structure 
formation for F co h = 0, we expect that AEi nt (and con- 
sequently /3) also scales linearly with F co h to first order. 
This implies 



Ao 



coh 



1/2 



(ii) 



where s is a constant with units of stress. 

This toy model makes a number of important simpli- 
fications. Foremost is the fact that near the aperture 
the 'real' streams are more complicated than we have 
assumed. Much of this results from the presence of ra- 
dial velocity gradients (both dv z /dr and dv r /dr), which 
change the initial kinetic energy and the initial rate at 
which energy is lost. 

Nonetheless, in comparing our idealized model with 
data from the simulations we find good agreement. In 
Fig. [9] we plot the droplet lengthscale measured at the 
aperture as a function of the cohesive force along with 
a fit to Eq. ITT] with s as a free parameter. As can be 
seen, the F c ' oh scaling of the model is consistent with 
behavior found in the simulations. We see that at ~ 100 
nN, the initial lengthscale is approximately equal to the 
particle size (200 /xm). Intuitively, one might expect this 
model to break down once the initial lengthscale drops 
below the particle size, and the fact that the force at 
which this occurs is very near to the crossover identified 
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FIG. 9: (Color online) Droplet lengthscale at the aperture. 
Data are for the same grain diameter and aperture size as in 
Fig. [2] Measured droplet separation at aperture (red circles) 
and toy model fit (solid line). The single fit parameter has 
the value s = (2.05 ± 0.03) x 10~ 2 nN/mm 2 . 

in Fig. [6] is highly suggestive. The fit corresponds to 
s = (2.05 ± 0.03) x 10" 2 nN/mm 2 . Referring back to 
Eqs. [lO] and [TT] we see that this quantity encapsulates 
the ratio of the initial kinetic energy Ki to the net energy 
lost AE int — W (divided by factors of F co h and Ao to 
give the appropriate units). Physically, it is a proxy for 
the stream's effectiveness to counteract its initial kinetic 
energy through cohesion and dissipation. While it would 
be preferable to have a larger range of data to fit to, this 
is not possible in our stream geometry. As can be seen 
in Fig. [4j the droplet forming regime gives way to the 
dripping regime quite rapidly just beyond the maximum 
force used in Fig. [9j and we are left with only about a 
decade in F co h and Ao before a new mechanism comes 
into play. 



II. CONCLUSIONS 

Our results show that freely falling granular streams 
provide a highly sensitive probe of minute, local-level in- 
teractions between granular particles. Gravity stretches 
and cools the granular stream, similar to what happens 
during the expansion phase in molecular jets. This makes 
it possible to distinguish between situations with and 
without attractive particle interactions, irrespective of 
the coefficient of restitution. Attractive interactions that 
are weak (nanoNewtons) , short-ranged (tens to hundreds 
of nanometers) and thus hard to quantify in-situ with 
other techniques are easily detected by tracking the for- 
mation of particle clusters or droplets. In particular, our 



findings demonstrate that the average number of parti- 
cles in a cluster or droplet, given by the saturation net- 
work size, is directly related to the cohesive strength. 
This also implies a simple experimental means for detect- 
ing small amounts of interparticle cohesion without the 
need to follow the stream in the co-moving frame: From 
images taken in the large- z limit the number of particles 
per droplet can be estimated, providing an upper bound 
on the cohesive energy and making these streams a po- 
tentially useful tool for other areas of granular physics 
where cohesive forces are known to play important roles 
but are difficult to quantify, such as industrial fluidized 
beds [7] and accretion in protoplanetary disks [30] . 

Nearly all aspects of stream behavior are affected by 
the free-fall velocity profile imparted on the stream at the 
aperture. This observation suggests that the breakup 
mechanism itself may be a consequence of this gradi- 
ent. With this in mind, we developed a toy model in 
which droplets are formed as the kinetic energy associ- 
ated with the gradient and the work done by neighboring 
droplets competes with the energy stored and bonds and 
lost to dissipation. This model predicts that adjacent 
droplets should separate with constant velocity, which is 
confirmed in the simulations. Using this model to deter- 
mine the length of droplets requires some knowledge of 
how the energy is transferred between the initial and final 
states. Predicting this is difficult because energy is lost 
to inelasticity, friction, and the rearrangement of bonds, 
all of which depend on the microscopic dynamics dur- 
ing the breakup process. Nonetheless, we expect that, to 
first order, these terms should just scale linearly with the 
cohesive force. This results in the initial droplet length- 

1 /2 

scale scaling like F c ' oh . This prediction is consistent with 
our data for droplet forming streams in the range of data 
available to our stream geometry. 

Our results indicate that in two very important ways 
granular breakup differs from inviscid Rayleigh-Plateau 
breakup. First is the fact that in the liquid case the 
lengthscale over which breakup can occur is independent 
of the strength of the interaction; surface tension is only 
required to start the process. For granular breakup, our 
results show that the size of droplets grows with the cohe- 
sive energy, indicating a competing mechanism is present. 
Second, in granular breakup dissipation plays a primary 
role, while in liquid breakup dissipation is not even nec- 
essary. In inviscid liquid breakup, droplets actually end 
up with more kinetic energy than they start with. This 
is a consequence of energy conservation; in order for a 
liquid stream to reduce its surface energy, its kinetic en- 
ergy must increase. In granular breakup, it is clear that 
fully formed droplets have virtually no kinetic energy in 
the co-moving frame. 

In a liquid stream, it is the thermal nature of the 
molecules which allows them to explore phase space 
and ultimately reduce the surface energy of the stream. 
Macroscopic granular systems, on the other hand, are 
athermal, and if we imagine a well-connected column of 
grains with no macroscopic velocity gradients, we would 
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never see breakup. In the presence of a velocity gradient, 
grains have the opportunity to rearrange in a systematic 
way that maintains high density, allows access to energet- 
ically favorable configurations (by increasing the number 
of bonds), and leads to energy loss via dissipation. In 
this sense, the kinetic energy from the velocity gradient 
in a granular stream plays a similar role to the thermal 
energy in a liquid stream in that it provides a way for the 
stream to explore otherwise unaccessible configurations. 
However, in the absence of external energy input this 
kinetic energy is quickly lost as the grains explore new 
configurations and is essentially zero once the drops sep- 
arate. Combined with the extremely low effective surface 
tension, this prevents the granular droplets from assum- 
ing the smooth spherical shape of liquid drops. 

While our model and predictions can account for many 
of the trends and behaviors seen in both simulations and 
experiments, they also raise a number of important ques- 
tions. One of these is how the break-up is changed, or 
if it even occurs, with other interactions such as conser- 
vative potentials or frictionless particles, since changing 
these will significantly affect how the particles can redis- 
tribute their kinetic energy. Also, while our model pro- 
duces scaling for F co h consistent with the simulations, it 
is not able to determine what happens between the initial 
and final states, nor is it able to explain why the aspect 
ratios of fully formed droplets are largely independent of 
the stream parameters. As a final question, one has to 
wonder how the initial, local particle packing inside the 
stream near the aperture affects the breakup since the 
collision dynamics and the stream's ability to create new 
bonds will be strongly affected by this parameter. 

Finally, the fact that already minute amounts of cohe- 
sion can give rise to cluster or droplet formation provides 
a new perspective for interpreting other recent experi- 
ments. This includes the liquid-like breakup of highly 
collimated jets that rise up vertically after a heavy object 
impacts fine, loose granular material [31 ! and the clus- 
tering inside thin granular sheets that form after a dense 
granular stream hits an obstacle [32]. As in the granu- 
lar streams, heterogeneities become visible while the jets 
or sheets expand after a brief stage of intense collisional 
interactions. Inelasticity was suspected as the source, 
but our findings show instead that the observed dynamic 
structure formation is indeed the vestige of residual at- 
tractive interactions. 



III. APPENDIX 

Here we describe the technical aspects of the simula- 
tions in more detail. Three-dimensional molecular dy- 
namics simulations of both hopper and stream are per- 
formed using Itasca PFC3D software. The system con- 
sists of a cylindrical hopper (diameter Dh) with a circu- 
lar aperture (diameter Do) at the center of its base. The 
hopper geometry is designed to produce constant flow 
with the smallest number of grains necessary [33] . Main- 



TABLE I: Simulation parameters 
Parameter Value 



d 


100 /im, 200 /im 


P 


2.5 kg m -3 


Do 


2.0 mm, 3.0 mm 


D H 


5.0 mm, 7.5 mm 


N hop 


7xl0 4 , 3xl0 4 


kn i kt 


500 N/m 




0.50 


lc 


100 nm 



taining D H > 2.5L> , D H > 30d + A), and a fill height 
H > 1.2Dh insures particle flow that is insensitive to 
the addition of further grains. To establish the proper 
fill height, the simulations maintain Nhop particles inside 
the hopper. Values for the simulation parameters used 
are given in the table below. For d,Do,Dn and Nhop 
the first value listed refers to the streams shown in Fig.JT] 
and the second to the images and data shown in Figs.[2][9] 
To initiate a simulation, grains are rained into the hop- 
per and allowed to settle under gravity into a granular 
bed of height « H before the aperture is opened (inci- 
dentally, we believe that the systematic deviations in the 
e = 0.61, 

F C oh — 1000 nN, 1100 nN streams may have 
their source in insufficient system initialization). Once 
the flow starts, grains are added to the top of the bed 
to keep Nhop and therefore H constant on average. The 
average particle velocities at the aperture produced by 
the simulations closely matched those found in experi- 
ments [18] with similar geometries (e.g., v z « lOcm/s for 
d =200 /im and D () =3 mm). 

The repulsive force between grains is modeled using 
a linear spring-dashpot with both normal and tangen- 
tial damping as well as static friction, following [25] . 
In this model, the normal contact force is given by 
F n = k n S n + 7nm e ff5 n , where S n = d - |r* - r^- 1 is 
the overlap between the grains i and j, k n the normal 
stiffness, 7 n is the damping constant for normal motion, 
m e ff = m/2 is the reduced mass, and m = ^npd 3 is the 
grain mass (p =2.5 kg/m 3 ). The tangential contact force 
F t has the same functional dependence, but with tangen- 
tial displacement 5 t , tangential stiffness k t and damping 
constant j t in place of 5 ni k n and j n . Friction is in- 
corporated by truncating F t in order to satisfy the yield 
criterion F t < pF n . 

Head-on collisional energy loss was parameterized by 
a coefficient of restitution e n = v' n /v n , where v n and v' n 
are the normal components of the colliding grains' rela- 
tive velocities immediately before and after the collision. 
This linear-dashpot mo del leads to a r estitution coeffi- 
cient e n = exp [— 7T7 n / ^ Ak n /m e f f — 7^], independent of 
impact velocity (note that we define e n as a 'bare' restitu- 
tion coefficient that does not take into account cohesion) . 
Similarly, a coefficient of restitution for the tangential 
motion e tl with k t and j t replacing k n and 7 n , captures 
losses due to sliding. For the simulations presented here, 
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we set k n = k t and j n — j t so that there is the same 
coefficient of restitution e n = e t = e for normal and tan- 
gential motion. The elasticity is varied by changing only 
the damping constants j n and j t while keeping fc n , k tl 
and m e ff constant. 

Cohesion is included by adding a force of constant mag- 
nitude F co h that turns on after two grains have collided 
and drops to zero outside a thin shell of thickness of 
l c around each grain. Physically, this model represents 
a force that turns off over a lengthscale that is much 
smaller than the particle size, e.g. van der Waals inter- 
actions or the rupture of a liquid bridge [26]. As a result, 
this is a hysteretic interaction and completely separat- 
ing two grains initially touching costs a cohesive energy 
W co h = J^ lc Fcohdr = F co hl c . We note that our model 
for cohesion between grains does not capture more sub- 
tle, higher order effects caused by, for example, a velocity 
dependent coefficient of restitution [see for example [TT]]. 

Droplet widths and heights are calculated for net- 
works with 30 or more particles in streams with F co h > 
100 nN. Effects of stream initialization (e.g. plug for- 
mation) are avoided by excluding the 1000 droplets far- 
thest from the aperture. To ensure droplets are fully 
formed/comopletely separated, only networks whose cen- 
ter of mass locations are greater than 20 cm from the 
aperture are included. The droplet height is defined as 
the maximum extent of a network along its principle axis. 
The droplet width is taken to be the mean radial extent 
from the principle axis of each droplet averaged over all 
angles (partitioned into 36° bins). 

In our system, the heterogeneous nature of the stream 
made it difficult to define a coarse grained velocity field. 
We instead define the mean squared velocity deviations 
((Av) 2 ) for each grain as the standard deviation of its ve- 



locity relative to the velocities all other neighbors within 
a small region, which we set at 1.5d. If a grain has less 
than three neighbors within this region it is excluded, 
thus ensuring that isolated grains do not contribute. The 
velocity deviations are calculated for each grain at time 
steps 0.2 ms apart, and then binned into 100 fim wide 
bins along the axial direction to make the plots in Fig. [5j 
To calculate the normalized average contact number 
and average network size, we define grains to be in con- 
tact if the separation between their surfaces was less 
than the thickness of the cohesive shell l c . The data in 
Fig.|6ja) are calculated for each grain for intervals 0.2 ms 
apart and then binned into 1 mm wide along z. Data in 
Fig. J6jb)-(c) were also calculated for each grain in inter- 
vals 0.2 ms apart but binned into 5 mm wide bins along 
z. 
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